A dynamic localization of 2D electrons at mesoscopic length scales 
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We have investigated the local magneto-transport in high-quality 2D electron systems at low 
carrier densities. The positive magneto-resistance in perpendicular magnetic field in the strongly 
insulating regime has been measured to evaluate the spatial concentration of localized states within 
a mesoscopic region of the samples. An independent measurement of the electron density within the 
same region shows an unexpected correspondence between the density of electrons in the metallic 
regime and that of the localized states in the insulating phase. We have argued that this correspon- 
dence manifests a rigid distribution of electrons at low densities. 

PACS numbers: 73.21.-b, 73.20.Qt 
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The nature of localization of electrons in low dimen- 
sional systems is a matter of continuing debate. The 
enhanced effect of Coulomb interaction in such systems 
gives rise to quantum many-body ground states which, 
depending on the extent of disorder, may display a wide 
variety of excitation spectra and real-space distribution 
of charge. From extensive theoretical |l| and experimen- 
tal [2| investigations, it is now generally agreed that at 
zero temperature (T = 0), a gap would appear in the 
single particle density-of-states (DOS) at the Fermi en- 
ergy (Ep), even though the form of this gap remains 
controversial. On the other hand, direct experimental 
determination of the spatial distribution of the localized 
states is significantly rare. Theoretically, several forms 
of the distribution have been predicted in the absence of 
disorder, which include Wigner crystals 0, bubbles or 
striped phases |4|. However, when disorder is present, 
strong localization is associated with the pinning of elec- 
trons to a static, randomly distributed impurity sites. 
In this paper, we have addressed this issue by directly 
probing the spatial distribution of localized states within 
a mesoscopic region of low-density 2D electron systems 
(2DES) at various levels of background disorder. We find 
that even in the strongly insulating regime, the areal den- 
sity of these states (n£) changes with the electron den- 
sity (n s ), and hence are "dynamic" in nature. When 
analyzed in the framework of percolation theory, a direct 
correspondence between nj, and n s could be established, 
implying an unexpected rigid distribution of the electrons 
that is independent of the topology of background disor- 
der, but depends only on Ep of the system. 

In disordered 2DES, the strongly localized regime is 
identified by a longitudinal resistivity p 3> h/e 2 , and an 
activated (nearest-neighbor) hopping or a variable-range 
hopping transport, where p(T) ~ po exp (Tq/T) p , Tq and 
p being the relevant energy scale and hopping exponent 
respectively. In a weak magnetic field (B), applied per- 
pendicular to the plane of the 2DES, the asymptotic be- 
havior of the wave-function of an isolated impurity state 
changes to ^(r) ~ exp (— r/£ — r 3 £/24A 4 ), where £ is 
the localization length and A = y/hfeB ^> £ 5J. The 



B-induced compression of the wave-function leads to a 
strong positive magneto-resistance (MR), which is ex- 
pressed analytically as p(T, B) = p(T, B = 0) exp (aB 2 ). 
An expression for the geometry-dependent factor a has 
been derived from the shift in percolation threshold, 
which in the limit of narrow bandwidth can be expressed 
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where At, is a model-dependent constant of order unity, 
whose precise value will be discussed later. Eq.Qprovides 
a mechanism of estimating til from the weak-field MR 
in strongly localized systems, and has been used in the 
context of hopping transport in the impurity band of 8- 
doped GaAs flfl. In the present work we have used this 
technique to calculate within a mesoscopic area in 
low-disorder 2DES, and compared it to the absolute value 
of n s within the same region. 

As shown in the schematic of Fig. la, mesoscopic seg- 
ments of high- mobility Si <5-doped GaAs/AlGaAs het- 
erostructures were created by the intersection of a narrow 
wet-etched mesa (width: W) and a transverse Au/NiCr 
metallic gate (width: L). The doping density ns = 
2.5 x 10 12 cm~ 2 was kept constant for all samples, while 
the background disorder was tuned by changing the sepa- 
ration (8 sp ) of the dopant layer and the GaAs-AlGaAs in- 
terface. The value of n s within the active region depends 
on the gate voltage (V g ) through the specific capacitance 
Cq. Within the simple capacitor model, Cq ~ eoe r /ed s , 
where d s is the depth of the 2DES from the surface. For 
this experiment we have chosen mesoscopic samples from 
two heterostructures, which differ in S sp and d s , but iden- 
tical in all other structural and geometrical aspects. The 
relevant parameters are given in Table ID Note that sam- 
ple A07 is more disordered than A78 due to a smaller 
S sp , and hence has lower electron mobility (p). All ex- 
periments were performed at T « 0.3 K unless specified 
otherwise. In order to minimize the effect of contact and 
series resistances, the electrical measurements were car- 
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TABLE I: Relevant parameters of the samples used. p(0) and 
n s (0) are as-grown mobility and electron density respectively. 



Sample 


M (0) 

cm 2 /V-s 


n 3 (0) 
cm -2 


5sp 

nm 


ds 
nm 


L x W 
/jmxjjm 


Co 

chi^V -1 


A07 


0.6 x 10° 


2.9 x 10 11 


20 


120 


2x8 


57 x 10 iu 


A78 


1.8 x 10 6 


2.1 x 10 11 


40 


290 


2x8 


24 x 10 10 



ried out only in 4-probe geometry with a low-frequency 
(~ 7.2 Hz) ac excitation current of ~ 0.01 — 0.1 nA. The 
maximum measured resistance was limited to ~ 1.5 Mf2 
to avoid any systematic error arising from stray capaci- 
tances. 

The first part of the experiment evaluates the 
V g — dependence of n s when n a is sufficiently high such 
that p <C h/e 2 and quantum Hall-based methods are 
applicable. However, due to the restricted dimensions 
of the active region, in particular the smallness of L, 
conventional Hall or Shubnikov-de Haas measurements 
are experimentally unfeasible in these samples. Hence, 
we have used an alternative technique based on the re- 
flection of edge-states, which was originally developed to 
study backscattering of edge-states by well-defined po- 
tential barriers The technique is schematically rep- 
resented in Fig. la. Briefly, when the perpendicular B 
corresponds to vq edge channels in the ungated part of 
the Hall-bar, a negative V g will cause — v of the chan- 
nels to get reflected at the gate, where v is the filling- 
factor within the active region. As a function of V g , this 
leads to plateaus in the four-probe resistance (i?i2,34) at 
-Ri2,34 = {h/e 2 ){\/v — \/vq), when the contacts (1-4) in 
Fig. la are close to ideal with no reflection and perfect 
transmission. This is illustrated in Fig. lb, where the 
V g — dependence of i?i2,34 is measured in sample A78 at 
a constant B corresponding to v$ — 10. Three plateaus 
corresponding to v = 8, 6 and 4 appear nearly exactly 
at the theoretically expected values (horizontal dashed 
lines). The V g at the center of these plateaus then cor- 
responds to an electron density of n s = eBv/h, while 
the width of the plateau represents the maximum uncer- 
tainty. For both samples, we have repeated this measure- 
ment for different values of Uq (within spin degeneracy) . 

Fig. lc shows n s as a function of V g for samples A07 
(empty symbols) and A78 (filled symbols). An advan- 
tage of this method is the negligible uncertainty in the 
y-direction since n s is derived from a known integer (y). 
At a sufficiently negative V g both samples display the 
expected linear dependence of n s on V g . In A07, the 
deviation from linearity at higher V g can be explained 
from screening of V g by accumulation of carriers in the 
dopant layer [9j . When the linear part is fitted with n s — 
C s {V g - V° s ), the slope C s was found to be » 55.5 x 10 10 
cm^V- 1 for A07 and « 23.0 x 10 10 cm^V- 1 for A78, 
agreeing satisfactorily with corresponding values of Co- 
The intercept V° s was found to be V° s = -0.917 ± 0.004 



V and -0.804 ± 0.002 V for A07 and A78 respectively. 

We now consider the strongly localized regime, where 
p ^ h/e 2 . As shown in the insets of Fig. 2, at B = 0, 
this regime corresponds to V g < —0.86 V for A07 and 
V g < —0.74 V for A78. Increasing B from zero initially 
results in some negative MR in both samples, followed by 
an exponential rise in p(B) at higher B. The dependence 
of In p on B 2 at different V g for both samples is shown 
in Fig. 2. In order to ensure consistency, the analysis is 
limited to the regime of V g for which p(B = 0) > 2h/e 2 . 
From Fig. 2, within the measured range of p, lap varies 
linearly with B 2 over a field range of 1.2 T < B < 3 T for 
A07 and 0.5 T < B < 1.3 T for A78. The lower limit of B 
is determined by the extent of negative MR, which arises 
from an interference among the hopping paths |l0j . The 
upper cut-off in B probably arises from the strong-field 
condition r <C A 2 /£ for the asymptotic form the wave- 
function ip(r) ||- This is not fully understood at present 
and will be discussed in a future publication. 

The most striking aspect of the data shown in Fig. 2 is 
the continuous decrease in the slope a of In p — B 2 traces 
as V g is made more positive (i.e. as n s is increased). This 
is observed in both samples (Fig. 3). Note that from 
the expression for a in conventional hopping framework 
(Eq.^l, a is expected to increase with increasing V g since 
the localization length £ increases rapidly as the system 
becomes more metallic (til remains constant since it is 
assumed to arise from frozen impurity distribution) . This 
has been confirmed in hopping MR in the Na + impurity 
band of Si-MOSFET's, where a was found to increase as 
the Fermi energy was swept towards band half-filling |ll| . 
In our case, we also note that the absolute magnitude of 
a is about factor of ~ 2 larger in A78 than A07 over the 
same range of p. 

If we assume the validity of Eq.^ the counter- intuitive 
trend of a in Fig. 3 forces ul to be V^-dependent. Fol- 
lowing Ref. pj, since p 3> h/e 2 , we first replace £ of 
the hydrogenic wave-function of the strongly localized 
states by a* B , where a* B is the effective Bohr radius of 
electrons. However, an evaluation of til from Eq. Q]also 
requires the magnitude of Al- If the system consists of 
only two isolated localized states, one obtains Al sim- 
ply from the overlap integral as Al = 1/12 k 0.083 [l2j |. 
However, in the hopping regime there are a large number 
of states, and the relation Al — (N c /n) 3 / 2 /12 can be ob- 
tained following Nguen's analysis of the i?-induced shift 
in percolation threshold. Here N c is the dimensionality- 
dependent number of bonds per site in the random resis- 
tor network 6]. iV c = 4.5 or 2.7, for zero or finite width 
of the impurity band with respect to ksT 0,0. This 
gives Al « 0.143 and 0.0664 respectively. Note that all 
the estimates of Al are well within the order of unity, 
and hence only has a fine tuning effect on the absolute 
magnitude of ul- In the subsequent analysis we have 
used Al = 0.0664, even though the effect of other values 
of Al has also been discussed. 
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Fig. 4a shows the V g — dependence of ul evaluated us- 
ing Eq. ^and the experimentally obtained a. Note that 
for both samples displays an approximately linear de- 
pendence on V g . A fit of the form n L = Cl{V 9 — V° L ) 
gives C L w 55.7 xlO 10 cm^V" 1 and 23xl0 10 cm^V" 1 , 
and V° L k -0.918 ± 0.003 V and -0.806 ± 0.005 V for 
A07 and A78 respectively. 

The calibration of n s with V g now allows us to identify 
an unexpected correspondence between ul and n s for 
both samples. Firstly, the slopes C s and Cl agree with 
each other as well as the expected specific capacitance 
Co. Secondly, the pinch-off voltages V^ L and Vg S are also 
equal within the experimental accuracy. To confirm this 
correspondence, we have plotted til and n s together as 
functions of AV g , where AV^ = V g — V° L for til and V g — 
Vg S for n s . Fig. 4b clearly shows ul to be a continuation 
of n s in the low-density, strongly- localized regime in both 
samples. While the best alignment could be obtained 
using Al = 0.0664, Fig. 4b also illustrates the position 
of iii relative to n s for A07 with Al — 0.143 (dotted 
line) and Al — 0.083 (dashed line). 

The correspondence between n l and n s clearly demon- 
strates a distribution of localized states that is indepen- 
dent of the topology of disorder, but depends only on 
the density of electrons. Such a possibility has been 
discussed extensively in the context of Wigner crystal- 
lization in low-density systems. Recent quantum Monte 
Carlo simulations have indicated the possibility of sta- 
bilizing a Wigner-like rigid distribution of electrons at 
relatively high electron densities in presence of optimum 
disorder [14|. Microwave absorption studies have also 
claimed possible observation of Wigner solid in the frac- 
tional quantum Hall regime at n s as high as ~ 6 x 10 10 
cm~ 2 |l5j. Transport in such systems has been discussed 
in terms of hopping of vacancies or interstitials, or the ex- 
citation of such species to the mobility edge |16l Il7j | . In 
this experiment, even though we have probed the spatial 
distribution of localized states in a more direct manner, 
we cannot, at least at present, make a definite statement 
on any order in the distribution. This would require fur- 
ther understanding of the geometrical parameters, in par- 
ticular Al, which is presently under investigation. 

The temperature dependence of a provides further in- 
sight to the nature of the "impurity" band and corre- 
sponding DOS (£>l). In the inset of Fig. 3b we have 
shown a typical T— dependence of a, corresponding to 
V g = —0.760 V for A78. Over the observed range of T, 
a is nearly constant. In the percolation framework, this 
confirms the validity of Eq^for our system and allows us 
to express the bandwidth W as W < fc^Tln (rj c ), where 
rj c determines the resistivity at the percolation threshold, 
p = Po exp (77c). To get an estimate of W, we use the ob- 
served p/po « 6 at V g = -0.76 V and T — 0.3 K. {p ~ 
h/e 2 was obtained from T— dependence of p.) This gives 
W < 0.05 meV, and from the observed til ~ 1 x 10 10 



cm~ 2 at the same V g , we get D L ~ n L /W w 2.1 x 10 14 
cm _2 eV _1 . Interestingly, this is nearly an order of mag- 
nitude larger than even the spin-degenerate free-electron 
DOS in GaAs, m* /irti 2 = 2.7 x 10 13 cm- 2 eV" 1 . Sim- 
ilar observation was reported in the context of trans- 
port near the mobility edge in the inversion layer of Si- 
MOSFETs [l8j. While a shift in the mobility edge as a 
function of Ep was suggested, in our case possibility of 
an interaction-induced enhancement of the effective mass 
(to*) cannot be ruled out. 

An important aspect of the solid phase is the spin 
of individual localized states. Rigid electronic distri- 
bution at til as large as 1 — 5 x 10 10 cm~ 2 will def- 
initely involve strong exchange interaction. Both fer- 
romagnetic and anti-ferromagnetic ground states have 
been predicted depending on the nature of ring exchange 
in this phase 0, 0|. Indeed, possible observation of 
a spontaneously spin-polarized phase, in particular at 
lower dimensions, is now being reported in several exper- 
imental investigation in the low-density regime |16l |20j . 

In conclusion, magneto-transport in the strongly lo- 
calized, low-density regime of high-quality 2D electron 
systems provide evidence of a dynamic, gate voltage- 
dependent distribution of localized states. Over a meso- 
scopic region of the 2D systems, the spatial density of the 
localized states was found to agree with the estimated 
density of electrons at the same gate voltage. This corre- 
spondence is argued to be manifestation of a many-body 
rigid distribution of electrons, driven by strong Coulomb 
interaction. 
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FIG. 1: (a) A schematic of the device structure. Contacts 1 and 2 are for current probes, while 3 and 4 used for voltage probes, 
(b) Illustration of gate voltage (Vg) dependence of the 4-probe resistance 7?i2,34 at a constant magnetic field (B = 0.742 T) 
corresponding to an integral filling factor uq (= 10) in sample A78. The dashed lines are the expected values of -Ri2,34 at 
different filling factors v within the gated region (see text), (c) n s — eBu/h, at values of V g 's corresponding to the center of 
plateaus in i?i2,34. The horizontal error bars represent the width of the plateaus in V g . The solid lines are fits to the linear 
part of the data. For both samples, different symbols correspond to different integral value of vq at which i?i2,34 — V g sweeps 
were recorded. 




FIG. 2: Dependence of resistivity p on B 2 for different values of gate voltage (V g ) recorded at T — 0.3 K. Inset: V g — dependence 
of p at zero and finite B. (a) Sample A07, and (b) Sample A78. 




FIG. 3: Gate voltage (V 9 ) dependence of the slope (a) of the linear region of In (p) — B 2 data shown in Fig. 2. The data was 
obtained at T = 0.3 K. The fit uncertainty increases, and hence the error bar on a, as p decreases to ~ h/e . (a) Sample A07, 
and (b) Sample A78. Inset of (b): the typical temperature dependence of a obtained at V g = —0.760 V for A78. 
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Fig.4 




FIG. 4: (a) The spatial density of localized states (til) as a function of V g calculated from Eq. 1 and the experimentally 
obtained a. Solid lines denote linear fit. (b) Correspondence between iil and n s - For comparison, til, n s ~ AVg are plotted 
in log-log scale, where AV g — V g — V g ° L for til (empty symbols) and V g — V ga for n s (filled symbols). For A07, the solid, dashed 
and dotted lines represent estimates of ul assuming Al = 0.0664, 0.083 and 0.143 respectively. 



